jsMath
Uvjetovanost Hilbertove matrice
def hilbert_matrix(n): return matrix([[1/i for i in [j..j+(n-1)]] for j in [1..n]])
def hilbert_cond(n): A = hilbert_matrix(n) return A.norm()*A.inverse().norm()
hilbert_cond(5)
476607.250243
hilbert_matrix(30).inverse().norm()
2.12827697168e+43